Axion Dark Matter and Cosmo-logical Parameters 
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We observe that photon cooling after big bang nucleosynthesis (BBN) but before recombination 
can remove the conflict between the observed and theoretically predicted value of the primordial 
abundance of 7 Li. Such cooling is ordinarily difficult to achieve. However, the recent realization that 
dark matter axions form a Bose-Einstein condensate (BEC) provides a possible mechanism, because 
the much colder axions may reach thermal contact with the photons. This proposal predicts a high 
effective number of neutrinos as measured by the cosmic microwave anisotropy spectrum. 

PACS numbers: 95.35.+d 



INTRODUCTION 

The agreement between observations and the BBN pre- 
dictions for the primordial abundances of light elements 
is often touted as a triumph of the standard ACDM cos- 
mological model. Under the assumption that there are 
three neutrino species, BBN as a theory requires essen- 
tially a single input: the baryon-to-photon ratio, usually 
given by the parameter r/io = 10 10 ris/n 7 [![. If one takes 
7710 to be 6.190 ± 0.145, in accordance with the latest 
Wilkinson Microwave Anisotropy Probe(WMAP) results 
the inferred primordial abundances of the majority 
of the light elements (D, 4 He, 3 He) are remarkably con- 
sistent with BBN predictions, save one exception: that 
of 7 Li is approximately two to three times less than what 
the theory predicts. The discrepancy is deemed statisti- 
cally significant, and there is so far no widely accepted 
explanation for the anomaly. In the literature, this is 
referred to as the "Lithium Problem" . 

One of the most difficult issues involved in testing BBN 
is how reliably to infer the primordial abundances of light 
elements from measurements that are available to us. 
Subsequent to BBN, the original relic abundances are 
all subject to further modification by complicated stellar 
processes. 7 Li, for example, can be both depleted and 
synthesized in stars, as well as produced by cosmic-ray 
nucleosynthesis. As such, the abundance of 7 Li is in- 
ferred primarily from absorption lines in the atmosphere 
of galactic halo stars with low metallicity, since these 
stars are very old and have experienced very little nu- 
clear processing (See [Hill for details). 

Although these post-BBN effects lead to considerable 
complication, they also open up many different avenues 
to explain the 7 Li anomaly. For many years, it has 
been hoped that better determination of nuclear param- 
eters will gradually narrow the discrepancy, though it 
was eventually realized that does not seem achievable |6j . 
Quite the contrary, it was found in [l| that improved data 
on the neutron life-time and the cross sections p(n,7)d 
and 3 He(a,7) 7 Be increases the predicted abundance of 
7 Li, worsening the disagreement. Revisions to stellar 
evolution, as a consequence of systematic errors in the 



effective temperature of the metal-poor stars 

0,1, 

and 

surface 7 Li depletion in the interior of stars due to some 
mixing or diffusive processes Q, have also been investi- 
gated as possible solutions, but are still considered con- 
troversial |l|. 

The fact that the nuclear reactions relevant to the pro- 
duction of both primordial and post-BBN 7 Li are now 
quite well understood has led to speculations that the 
anomaly might instead be caused by new physics. Many 
explanations have been proposed, such as the variation in 
time of the deuteron binding energy and of fundamental 
couplings [l(| 11 1, and the decay of a relatively long- 
lived particle in the context of supersymmetry [12| . At 
this point, none of these explanations have won general 
acceptance in the cosmology community. 

In this paper, we propose a way to remove the con- 
flict between data and theory for the abundance of Li, 
based on the cooling of photons between the end of BBN 
and decoupling. Processes that do this are difficult to 
come by. Indeed, typical processes arising from new 
physics tend to heat up the photons, modifying 7710 in 
the wrong direction [l(. However, the recent realiza- 
tion that dark matter axions form a BEC at approxi- 
mately 500 eV photon temperature [HI provides a pos- 
sible mechanism [3]. Essentially, the high occupation 
of axion modes with very low momenta greatly enhances 
the strength of their gravitational interactions, such that 
an exchange of energy between the photons and the much 
colder axions becomes possible. Photon cooling implies 
that 7710, bbn < ?7io,wmap, which has the effect of reduc- 
ing the production of 7 Li [5] . If thermal equilibrium be- 
tween the photons and axions is achieved, the 7 Li abun- 
dance is reduced by approximately a factor 2 (see be- 
low), alleviating the discrepancy and perhaps removing 
it altogether. However, our proposal predicts a higher 
abundance for D than present observations indicate and 
predicts that the effective number of thermally excited 
neutrino degrees of freedom is high: N c ^ = 6.77. 

Photon cooling by kinetic mixing with hidden photons 
was proposed in ref. fisl ]. 
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DARK MATTER AXIONS 

The axion was originally postulated to explain the ab- 
sence of CP violation in the strong interactions (l6| . It 
was later realized that the population of axions produced 
by the turn on of the axion mass during the QCD phase 
transition possesses the right properties to be cold dark 
matter (CDM) [T3] - First, they have the measured CDM 
density if the axion mass m is of order 10 _5 eV/c 2 [l8j | . 
Second, the axions thus produced are very cold: their 
average momentum is only of order the Hubble expan- 
sion rate (3 x 10 _9 eV/7i) when the axion mass effectively 
turns on and has been redshifting ever since to a mere 
10~ 17 mc today. Third, axions in this mass range interact 
very weakly through all forces other than gravity. These 
properties make axions one of the leading candidates for 
CDM. The other main contenders are weakly interacting 
massive particles (WIMPs) and sterile neutrinos. Hence- 
forth, we set h = c = 1 . 

Observationally it seems difficult to distinguish among 
the CDM candidates. However, it was recently real- 
ized that axions form a BEC through gravitational self- 
interactions at approximately a photon temperature of 
500 eV IH 14 | . The relaxation rate of cold dark matter 
axions through gravitational self-interactions is of order 

T a - iTTGnm 2 ? 2 (1) 

where n is the number density of cold axions and I ~ 
ti-^rr\ their correlation length. t\ ~ 2 • 10~ 7 s is the 
time at which the axion mass effectively turns on during 
the QCD phase transition, and a(t) is the scale factor. 
Eq. ([TJ is appropriate when the energy dispersion of the 
particles is less than their relaxation rate. We refer to this 
case as 'the condensed regime', to distinguish it from the 
more commonly encountered 'particle kinetic regime' de- 
fined by the condition that the energy dispersion of the 
particles is large compared to the relaxation rate. Cold 
dark matter axions are in the condensed regime after t%. 
The ratio T a (t)/ H(t), where H(t) is the Hubble rate, is 
of order 5 • 10 -7 at t\ but increases with time as a(i) _1 i, 
and reaches one at approximately 500 eV photon tem- 
perature. At that time the axions thermalize and form 
a BEC. Almost all axions go to the lowest energy state 
available. The correlation length grows and becomes of 
order the horizon. In the linear regime of evolution of 
density perturbations and within the horizon, the lowest 
energy state is time independent and no rethermalization 
is necessary for the axions to remain in the lowest energy 
state. In that case, axion BEC and ordinary CDM are in- 
distinguishable on all scales of observational interest [l3| ■ 
However, beyond first order perturbation theory and/or 
upon entering the horizon, the axions rethermalize to try 
and remain in the lowest energy available state. Axion 
BEC behaves differently from CDM then and the result- 
ing differences are observable. 



The study of the catastrophe structure of the inner 
caustics of galactic halos provides evidence that the dark 
matter is an axion BEC 13, 19[ . Briefly, there is a di- 
chotomy in the classification of the inner caustics in terms 
of their catastrophe structure, depending on the angu- 
lar momentum distribution of the infalling particles [20j . 
Axions in a BEC are in a state of net overall rotation 
and produce caustic rings, whereas ordinary CDM has 
an irrotational velocity field and produces tent-like caus- 
tics. There are several pieces of evidence for the existence 
of caustic rings at the predicted radii in various galaxies 
21 , 22 1 . It is shown in ref . [l9[ that the phase space struc- 



ture implied by the evidence for caustic rings is precisely 
and in all respects that predicted by the assumption that 
dark matter is a rethermalizing axion BEC. 

PHOTON COOLING 

The Lithium Problem refers to the mismatch between 
the observed and predicted abundance of primordial 7 Li 
by a factor 2 or 3 [J 0, Hi- For ^ 2.7, the predicted 
7 Li abundance increases with 7710- Hence, a cooling of the 
photons between the end of BBN and decoupling reduces 
the discrepancy. 

The gravitational fields of the cold axion fluid cause 
transitions between momentum states of other particle 
species present. For particles which are bosons or non 
degenerate fermions, the relaxation rate through gravi 



tational interactions with the cold axions is of order 14 1 



Ap 



(2) 



where uj is the typical energy of the particles and Ap 
their momentum dispersion. Eq. ([2]) generalizes Eq. (fTJ) 
to other species that are in the presence of the cold ax- 
ions. [Eq. (TT]) follows from Eq. ^ by setting uj = m 
and Ap = l~ x as is appropriate for the cold axions them- 
selves.] For photons to cool substantially it is necessary 
that energy is transferred from the photons to the low 
momentum highly occupied axion states and from those 
to the relativistic axion states. For both relativistic ax- 
ion states and for photons, Ap ~ lo and hence their re- 
laxation rate T r through gravitational interactions with 
cold axions is of order AnGnml. Using the Friedmann 
equation, one finds that T r /H oc a(t) before equality be- 
tween matter and radiation and remains constant after 
that. At equality, T r /H\ toci ~ ^(i C q)Acq- If i/t is order 
one at equality, the photons reach thermal equilibrium 
with the axions and hence cool. 

Gravitational interactions conserve particle number 
and therefore produce only kinetic (as opposed to chemi- 
cal) equilibrium between the species involved. Also, after 
500 eV photon temperature, the coupling between pho- 
tons and baryons is in the kinetic, rather than chemi- 
cal, equilibrium regime [231] ■ Upon cooling, the photons 



3 



that cannot be accommodated in thermally excited states 
enter the ground state, a plasma oscillation with zero 
wavevector. Since the photon chemical potential remains 
zero, the final photon spectrum is Planckian, consistent 
with observation. 

Eq. ^ does not apply to degenerate fermions be- 
cause of Pauli blocking. The cosmic neutrinos are semi- 
degenerate since they have a thermal distribution with 
zero chemical potential. Their thermalization rate is less 
than that T r of relativistic bosons. Since T r /H oc tnl oc 
t 2 a~ 3 (t), that ratio does not grow after equality. Since 
the relativistic axions may only reach thermal contact 
with the cold axions at equality and the neutrinos are de- 
layed relative to the relativistic axions, we believe it most 
likely that neutrinos remain decoupled from the axions, 
photons and baryons at all times. 

It is straightforward to determine how much the pho- 
tons cool if they reach thermal equilibrium with the ax- 
ions. Energy conservation implies pi n — p/, 7 + pf, a be- 
cause the contributions to the energy density of the initial 
axions and of the baryons are negligible. The ratio be- 
tween the final and initial photon temperature is thus 
(2/3) 1 / 4 . Since their number density is proportional to 
T 3 , we find: 

/2\ 3/4 

?7io,bbn=(-J r/i ,wMAP = 4.57 ± 0.11 (3) 

using 7710, wmap — 6.190 ± 0.145 [2]. Because the 7 Li 
abundance is proportional to rj 2 Q BBN in the range of in- 
terest, it is reduced by approximately the factor ~ 

0. 55. 

A number of authors proposed earlier that the dark 
matter is a BEC [HI] ■ The photon cooling described here 
may occur in those cases as well. If the particles are 
in the condensed regime, the relaxation rate is given by 
Eq. (fTJ), calculated with appropriate values for n, m, and 

1. However, in many of these proposals, the cosmological 
history of the dark matter particle is not known, render- 
ing the computation of the relaxation rate difficult. 

COSMOLOGICAL CONSEQUENCES 
Effect on light element primordial abundances 

Whether photon cooling by axion BEC solves the 
Lithium Problem remains to be seen. The data have been 
time dependent in addition to the usual uncertainties. In 
Fig. 1, we plot the value of 7710, bbn in the standard cos- 
mological model, labeled 'WIMP', and in the scenario 
described here, labeled 'axion', along with the values in- 
ferred from the observed light element abundances ac- 
cording to the review by G. Steigman in 2005 4], the re- 
view by F. Iocco et al. in 2008 [25| and a private commu- 
nication from G. Steigman updating his 2005 estimates 
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FIG. 1: Values of 7710, bbn inferred from the abundances of 
7 Li, D, 3 He and 4 He, and the predicted values in the standard 
cosmological model (WIMP) and in our proposal (axion) . The 
data inferred values are taken from refs. [2^] and [2^ |. 
The error bars indicate the 7710, bbn values consistent with the 
estimated 1-<t uncertainties in the observations. 



in the light of recent observations |26| . The error bars in- 
dicate the range of 7710, bbn consistent with the estimated 
1-cr uncertainties in the observations. The axion predic- 
tion agrees very well with the 7 Li abundance at the time 
of Steigman's 2005 review (?7 10 / L i = 4.50 ± 0.30). How- 
ever more recent observations indicate a lower primordial 
7 Li abundance, worsening the Lithium Problem. 

Perhaps more problematic is that a smaller 7710, bbn 
predicts an overproduction of D. Traditionally, D has 
been the prime choice as a baryometer among the light 
elements, due to its sensitivity to 7710, bbn and simple 
post-BBN evolution (abundance monotonically decreas- 
ing). The major drawback with D is that its abundance 
is inferred from a very small set of (seven) spectra of 
QSO absorption line systems [27] . Worse yet, these few 
measurements have a large dispersion, and do not seem 
to correlate with metallicity, obscuring the expected deu- 
terium plateau. Due to the various inadequacies in the 
D measurements mentioned, we have reservations about 
the common practice of attaching most significance on D 
in the comparison between data and BBN predictions. In 
comparison, 7 Li is inferred from a large number of mea- 
surements, which are more-or-less consistent. Also, since 
D is more easily destructible than 7 Li, it is conceivable 
that unknown stellar processes further deplete D. 

Finally the 3 He and 4 He inferred 7710, bbn values have 
large error bars and hence carry less statistical weight. 
The 4 He inferred value has increased recently (5.5 < 
7710. 4 Hc < 11 according to ref. [25| and 7.5 < ?7io. 4 Hc < 20 
according to ref. [26|) compared to its accepted value a 
few years ago, creating additional uncertainty. 
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Effective number of neutrino species 

After the axions are heated up and reach the same 
temperature as the photons, most of them are still in 
the ground state. The axions in the ground state behave 
as cold dark matter. The axions in the excited states 
contribute one bosonic degree of freedom to radiation. 
The radiation content of the universe is commonly given 
in terms of the effective number N e g of thermally excited 
neutrino degrees of freedom, defined by 



[2] E. Komatsu et al, Ap.J. Suppl. 192:18 (2011). 



Aad = p 7 [l +iV e ff^ f Jj- J ] 



(4) 



where p ra d is the total energy density in radiation and 
p~j is the energy density in photons only. The standard 
cosmological model with ordinary cold dark matter pre- 
dicts N c g — 3.046, slightly larger than 3 because the 
three neutrinos heat up a little during e + e~ annihila- 
tion. Taking account of the fact that not only is there 
an extra species of radiation (thermally excited axions) 
but also the contribution of the three ordinary neutrinos 
is boosted because the photons have been cooled relative 
to them, the proposed scenario predicts 



Prad 



P-j + Pa+ Pv 



Pi 



1 7 / 4 \ 3 3 

- + 3.046 - ( — ) - 



(5) 



which yields N cS = 6.77. 

At present, the measured values are smaller than this 
prediction. The WMAP collaboration found N c g — 
4.34+°;^ (68% CL) based on their 7 year data combined 
with independent data on large scale structure and the 
Hubble constant Q. An analysis 28 1 using the Sloan 
Digital Sky Survey (SDSS) data release 7 halo power 
spectrum found N e & = 4.8 ± 2.0 (95% CL). The At- 
acama Cosmology Telescope (ACT) collaboration finds 
29] N cS = 5.3 ± 1.3 (68%CL) using only their CMB 
anisotropy data and N e s = 4.56 ± 0.75 (68% CL) when 
combining that data with large scale structure data. The 
tendency for the measured values to be larger than 3.046 
has been taken sufficiently seriously to prompt propos- 
als for new physics involving extra neutrino species or a 
neutrino asymmetry (30l ) . The Planck mission is expected 
to measure N e g with much greater precision 31|. In so 
doing, it may shed light on the nature of dark matter. 

We thank Gary Steigman for extended discussions, and 
our referees for numerous suggestions for improvement. 
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